Quantitative Probe of Pairing Correlations in a Cold Fermionic Atom Gas 
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A quantitative measure of the pairing correlations present in a cold gas of fermionic atoms can be 
obtained by studying the dependence of RF spectra on hyperfine state populations. This proposal 
follows from a sum rule that relates the total interaction energy of the gas to RF spectrum line 
positions. We argue that this indicator of pairing correlations provides information comparable to 
that available from the spin-susceptibility and NMR measurements common in condensed-matter 
systems. 
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Introduction — The realization of degenerate atomic 
Fermi gases 0, 0, IE & IE IE nas opened new oppor- 
tunities for experimental discovery. The focus to date 
has mainly been on efforts to observe the condensation 
of atomic Cooper pairs to form a superfluid state simi- 
lar to the BCS state of electrons in a superconductor 0. 
Strategies for achieving observable pairing effects have so 
far hinged on the occurrence of strong attractive atom- 
atom interactions near a Feshbach resonance [E 0, ^| , 
in which a molecular bound state of one atom-atom scat- 
tering channel is close to the two-atom continuum thresh- 
old of another. The proximity of a Feshbach resonance 
can be adjusted by tuning a magnetic bias field, dras- 
tically altering the scattering behavior of atoms and al- 
lowing the s-wave scattering length to be varied over val- 
ues corresponding to effective interactions that are weak 
or strong, and repulsive or attractive. The scattering 
length, henceforth denoted by a , completely character- 
izes the interaction properties of atoms at low tempera- 
tures and low densities. The Feshbach resonance makes 

it possible to study one of the paradigms munii of 

fermion pairing theory, the BCS-BEC crossover, experi- 
mentally. 

Experimental groups have already observed the forma- 
tion of ther mal g ases 0, 0] and Bose-Einstein conden- 
sates |l7], [l8|, [l9j of diatomic molecules. Condensation of 
fermionic atom pairs on the attractive interaction side of 
the resonance, where there is no two-atom bound state 
and the analogy to the BCS-BEC crossover problem is 
closer, has also been reported [lEE]]- Because the BCS 
transition does not manifest itself strongly in the ex- 
panded density profile of the gas there is a need for quan- 
titative and direct measurements of pairing correlations, 
one that has motivated a large number of proposals. Sev- 
eral works focused on the change of light scattering due 
to the transition from the normal to the superfluid phase 
as a detection method [2E 0> HE Ell • Later, resonant 
laser light was proposed to induce tunnelling between the 
superfluid and normal state of the gas |2j| • The exper- 
imental realization of Bragg spectroscopy in a Bose gas 
llEl has inspired theoretical work on pairing effects in the 



dynamic structure function |29l |30|. Other interesting 
proposals include ones based on pairing induced changes 
in collective- mode frequencies Eil ll^lslil^lsl^lslij l. ro- 
tational properties of the gas [371 l38j | , the expansion of 
the gas|3J|, and atomic density noise-correlation prop- 
erties |40j. 

In this Letter we propose a more direct probe of pairing 
correlations that is similar to the spin magnetic suscepti- 
bility and nuclear-spin relaxation probes commonly used 
to detect electron pairing in condensed matter systems, 
and is able to detect pairing even when it does not lead to 
long-range coherence. We suggest a measurement of the 
cost in interaction energy when the number of Cooper 
pairs in the system is reduced by making the hyperfine 
state populations unequal. As we show below, this en- 
ergy change can be extracted from data obtained using 
the RF spectroscopy techniques that have already been 
developed by several experimental groups 0, 0, Eif . 
To illustrate the direct relationship between the hyper- 
fine population dependence of the interaction energy and 
pairing correlations, we compare the predictions of BCS 
theory for this quantity with its predictions for the more 
familiar magnetic susceptibility probe, which measures 
instead the dependence of total (interaction plus kinetic) 
energy on the same variable. 

Interaction Energy Sum Rule — We consider a gas of 
fermionic atoms that consists of a mixture of two hy- 
perfine states denoted by | |) and | J.). We assume the 
temperatures to be low enough so that only s-channel 
interactions, forbidden between atoms in the same hy- 
perfine state by the Pauli principle, are significant. The 
Hamiltonian of the gas is therefore H = Ho + Hi n t, where 
the noninteracting part is 
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and ip a is the fermionic annihilation operator for hyper- 
fine state The internal Zeeman energy of a hyper- 
fine state is denoted by e a , and for simplicity we have 
neglected any inhomogeneity of the magnetic field. In 
particular, this implies that we neglect the effects of the 
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magnetic trapping potential. We take the interaction be- 
tween unlike hyperfine states to be a contact interaction 
with strength V^, which should be chosen to produce the 
correct two-body scattering amplitude Q. With these 
assumptions, the interaction part of the Hamiltonian is 



H in t = V n /dx^j(x)^](x)^( X )V> T (x) 



(2) 



In the RF experiments one of the system hyperfine 
species (say J,)) is coupled to a spectator hyperfine state 
(|s)), and the number of atoms in the spectator state, N s , 
is detected as a function of the frequency of the RF-field. 
In the linear response limit N s is proportional to the 
rate of | J,) — > |s), transitions which we denote by 
We define the position of the associated RF spectrum 
absorption line as 
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(3) 



Using a formal golden-rule expression, I(uS) can be ex- 
pressed in terms of a two-particle correlation function of 
Fermion fields. It then follows from the fermion analog 

of sum rules derived in Refs. mill that 



(4) 



where V s -\ denotes the strength of the contact interac- 
tion between the spectator and 1 1) hyperfine states, and 
n a denotes the average density in spin state The 
shift of TuJ) from the bare line position huj Q = — e s dif- 
fers from the interaction energy per volume by a factor, 
^iT n t/(^it — ^st); wm ch can be held constant through 
the experiments and if necessary can be accurately de- 
termined by separate measurements. We conclude that 
the RF spectra enable a direct measurement of the inter- 
action energy density 
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Pairing and Interaction Energy in BCS Theory — The 
inverse spin-susceptibility of an unpolarized system of 
spin-1/2 particles may be expressed in terms of the de- 
pendence of free energy on spin-polarization: 



d8n 2 
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(6) 

where / to t is the total free energy per unit volume of 
the gas, n is the total density, and Sn = n-\ — is the 
spin density. It is well-known that the % s is strongly 
suppressed when atoms can gain energy by pairing. (\s 
vanishes as T — > in the BCS state.) For attractive 
atom-atom interactions, the energy cost of finite-spin po- 
larization has positive contributions from both interac- 
tion and kinetic energy. In the following we compare the 
spin-susceptibility with an alternate quantity that is de- 
fined in terms of the interaction energy alone and can be 
extracted from RF spectroscopy experiments performed 
for a series of hyperfine-state populations: 



-i _ d 2 e in t(n T ,nx) 
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As we show below, this quantity and the inverse spin- 
susceptibility provide similar probes of pairing correla- 
tions. 

We evaluate this interaction susceptibility using BCS 
theory from which it follows that 



eint(«T,^) 



m 



(8) 



where the dependence of the gap A = Vj-j, (V'i.V't) on 
temperature and hyperfine densities can be determined 
by solving the self-consistent mean-field equations. The 
mean-field Hamiltonian is 221 



f,[,i ls f ft 2 V 2 Anah 2 ni \ M , , t , . ( H 2 V 2 Anah 2 ^ \ , . . 

+ A^ (x) ^t (x) + AV;(x) ^ (x) _^_!!^^} , (9) 

1 * V n ml 



Note that the renormalization Vj-j, — » Anah 2 jm can be 
made at this stage in the Hartree mean-field potential. 
The chemical potentials of the two hyperfine states are 
denoted by jj, a , and are not necessarily equal, thus al- 
lowing for a density difference between the two hyperfine 



states. 

The partial densities are given by 

/dk 
^3 {\u k \ 2 N(fuv^ a ) + M 2 [1 - N(Hw ki - a )]} , 

(10) 
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where and are the Bogoliubov coherence factors, 
N(x) = \eP x + l] is the Fermi distribution function, 



and the | f ) quasiparticle dispersion is given by 



fruJ 



k.T 



(11) 



An identical expression, with the hyperfine labels inter- 
changed, applies for Swk.i- The Hartree-Fock mean- field 
shift is absorbed in the chemical potential via /i' Q = 
fi a — 4irah 2 n_ a /m. These equations for the densities 
need to be solved together with the BCS gap equation 

(2^) 3 2 v /[ek-(MT+Mi)/2] 2 + |A|2 V n ' 

Eq. 112|) contains a ultraviolet divergence that is renor- 
malized by using introducing the s-wave scattering length 
(a) between | f ) and | J.) hyperfine states by means of the 
Lippmann-Schwinger equation 
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dk 1 
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where we have introduced an ultraviolet cutoff fcA- 
follows H3 from Eq. O that 



It 



1 



n 



l 



2ak\ 



(14) 



For fcA — (lOOao)" 1 , and a = — 2000ao, we find that 
V\i — 0.07 x (Airah 2 jra). Note that although the spin 
densities and the BCS gap parameter are independent of 
the short-range properties of the interatomic potential, 
the interaction energy in Eq. (JSJ is not ■ 

Above the critical temperature, Tbcs — 
0.6T F e -^/2fep|o| ; A ->• so that 



Anah 2 



m 



(15) 



and xi^t s = 27r|a|ft 2 /m is temperature independent. For 
T < Tbcs the interaction energy is given by Eq. JSJ. 
Linearization of the gap equation implies that for T f 



Xint,s 



0.07fc F |a| 



— +27T\a\h 2 /m. (16) 



In Fig. n we plot Xi^± s anc l the inverse susceptibility 
vs. temperature, for the case of nj = = n/2, scat- 
tering length a = — 2000ao where ao is the Bohr ra- 
dius, and density n = 10 12 cm -3 . For these parameters 
Tbcs — 0.005 x Tp. The inverse interaction susceptibil- 
ity is greatly enhanced for temperatures below the BCS 
transition temperature, thus providing a clear signature 
of atomic Cooper pairs. 

Discussion and Conclusions — Although we have so 
far assumed a homogeneous Fermi gas, we believe that 
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FIG. 1: Plot of XhA s (solid line) and Xs -1 (dashed line) vs. 
temperature. Both susceptibilities are plotted in units of the 
quasiparticle density-of-states (dn/dfi). Note that the mean- 
field-theory interaction susceptibility is discontinuous at the 
critical temperature. 



our results also apply to the case of a trapped, and there- 
fore, inhomogeneous Fermi gas. There arc two main 
reasons for this. First, since the inhomogeneity in the 
density profile leads to broadening of the RF spectrum, 
the sum rule in Eq. @ is valid with the density n-f re- 
placed by a mean density 0> 03 ■ Second, because the 
inverse Fermi wave vector is much smaller than the har- 
monic oscillator length for the experimental systems of 
interest, the system may be treated within a local den- 
sity approximation j2^. Since the BCS gap parameter 
A will be largest in the center of the trap, the result in 
Eq. I|16fl should be evaluated at the center of the trap. 
With these two modifications, the results presented in 
this Letter should carry over, in a more than qualitative 
manner, to the inhomogeneous case. 

For strong attractive interactions, fc^H ~ 1, mean- 
field theory is not expected to be accurate. In particular, 
the superfluid transition temperature is expected to be 
limited by the loss of long range coherence rather than by 
the thermodynamics of pair formation. The thermody- 
namic probe we discuss here is sensitive to the occurrence 
of pair correlations and not particularly sensitive to the 
establishment of long range coherence. It should there- 
fore be able to detect the gradual development of pairing 
correlations with increasing interaction strength as the 
superfluid state is approached. We note that Bourdel et 
al. |47j | have measured the ratio of the interaction energy 
and kinetic energy of a Fermi gas by comparing expanded 
density profiles of an interacting gas of atoms with ex- 
pansion profiles of a gas at zero scattering length. In 
the weak-coupling limit such a measurement would pro- 
vide direct information on the temperature dependence 
of the interaction energy, since the kinetic energy is al- 
most independent of temperature in this case. In the 
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strong-coupling limit it is, however, not clear how the 
kinetic energy depends on the density difference, and it 
is, therefore, not obvious that a measurement of the ra- 
tio of interaction and kinetic energy for different hyper- 
fine state populations would provide a sensitive probe 
of pairing correlations in the gas in this limit. Finally 
we remark on the possibility of realizing inhomogeneous 
pair-condensate states in cold atom systems with un- 
balanced hyperfme state populations. These states could 
be detected by bringing the system to equilibrium in a ro- 
tating reference frame and visualizing their unusual |49| 
vortex-lattice structures. 
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